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I.  INTRODUCTION 


The  purpose  of  the  research  described  in  this  dissertation  is  to 
contribute  to  the  understanding  of  sound  attenuation  in  ocean  sediments. 
Although  the  measurements  were  performed  in  the  laboratory  on  laboratory 
(i.e.  artificial)  sediments,  they  will  aid  in  the  interpretation  of  in- 
situ  experiments  as  they  were  performed  under  controlled  conditions.  The 
specific  goals  of  the  experiments  presented  here  were  to  increase  under¬ 
standing  of  the  theories  of  sound  attenuation  in  sediments,  and  to  expand 
the  data  base  of  measurements  made  at  low  frequencies. 

Sound  waves  are  attenuated  far  less  than  light,  or  radio  waves  in  the 
1-2 

ocean,  and  have  been  used  extensively  for  communication,  navigation, 
and  remote  sensing.  In  deep  water,  the  transmission  loss  for  sound  waves 
is  well  understood,  as  the  ocean  bottom  plays  no  part  in  the  nroDagation 
of  acoustic  signals.  In  shallow  water,  the  transmission  loss  is  domi¬ 
nated  by  sound  absorption  in  the  ocean  floor.  This  transmission  loss  is 
the  key  to  understanding  the  prooagation  of  sound  waves  in  shallow  water. 

When  the  frequency  of  the  propagating  waves  is  high,  measurement  of 
the  attenuation  of  the  waves  is  not  difficult,  in  the  ocean,  or  in  the 
laboratory.  However,  at  low  frequencies,  the  long  wavelength  of  comores- 
sional  waves  causes  experimental  difficulties.  In  the  ocean,  the 
sediment  is  inaccesible,  the  propagation  path  is  inhomogeneous,  and 
information  on  layering  usually  has  to  be  derived  from  the  same 
experiments  as  the  attenuation  data.  In  the  laboratory,  a  large  scale 


2 


apparatus  is  required,  and  care  has  to  be  exercised  in  sample  prepara¬ 
tion.  The  propagation  velocity  of  shear  waves  is  considerably  smaller 
than  the  propagation  velocity  of  compressional  waves,  and  as  a 
consequence,  laboratory  and  in-situ  shear-wave  experiments  are  easier  to 
perform  than  compressional  wave  experiments.  This  dissertation  is  mainlv 

concerned  with  the  propagation  of  compressional  waves.  Some 

3 

investigators  feel  that  shear  wave  measurements  are  sufficient  to  under¬ 
stand  sound  propagation  in  sediments,  but  the  relationships  between  the 
propagation  parameters  of  shear  and  compressional  waves  are  based  udoo 
models  which  have  not  been  extensively  verified. 

This  dissertation  describes  a  method  for  measuring  the  attenuation  of 
compressional  waves,  at  low  frequencies,  in  laboratory  sediments.  The 
measurements  presented  here  were  performed  using  the  tube  method;  this 
method  has  been  widely  used  to  determine  sound  attenuation  in  other 
media.  The  tube  method  is  not  without  disadvantages,  the  effects  of  the 
tube  wall  on  the  propagating  wave  have  to  be  taken  into  account,  but  the 
method  does  eliminate  spreading  losses  and  it  is  relatively  easv  to 
generate  plane  waves  in  a  cylindrical  tube. 

Chapter  II  of  this  dissertation  presents  a  review  of  the  two  most 

widely  used  models  of  sound  attenuation  in  marine  sediments.  The 

4-1 1 

simplest  is  an  empirical  model  proposed  by  Hamilton,  the  other  model, 

to  which  most  of  the  chapter  is  devoted,  is  a  phenomenological  model 

12-18 

which  was  proposed  bv  Biot  as  a  general  model  of  sound  propagation 

3  19-21 

in  poro-elastic  media,  and  which  was  adapted  bv  Stoll  ’  for  use  on 


-A*-.- 


sediment  s . 
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In  chapter  III  the  physical  properties  of  the  sediments  used  in  the 
experimental  measurements  are  described.  The  calculations  of  thermal 
conductivity  and  specific  heat  for  sediments,  which  are  necessary  for  the 
calculation  of  the  attenuation  of  sound  at  the  tube  walls,  are  also 
described . 

Chapter  IV  is  devoted  to  a  description  of  the  experimental  system 


used  to  make 

the  sound  attenuation  measurements. 

Chapter 

V  contains  the 

analysis  of  the 

experimental 

data. 

The 

results  of 

the  experiments 

are  compared  with 

the  predicti 

ons  of 

the 

Hamilton  and 

Biot  models  of 

sound  attenuation 

in  sediments 

which 

are 

described  in 

chapter  II.  The 

final  chapter  summarizes  conclusions. 

'.-.I 
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The  functions  ber  and  bei,  which  derive  their  names  from  Bessel  real 
and  Bessel  imaginary,  are  zero  order  Kelvin  functions;  thev  are  related 
to  Bessel  functions  through^ 


1/2  B/2 

ber (W)+i  bei(W)=I  (i  W)=J  (i  '  W), 
o  o 


where  W  is  real. 


When  the  argument,  W,  of  of  the  viscosity  function  is  much  less  than 
one,  that  is  at  low  frequencies,  the  fluid  flow  in  the  pores  is 
Poiseuille  flow  and  the  viscosity  function  may  be  set  equal  to  one.  The 
viscosity  function  may  be  neglected  in  the  low-frequency  regime  of  the 
Biot  model.  The  changeover  from  low-frequency  behavior  to  high-frequency 
behavior  occurs  when  the  viscous  skin  depth  is  approximately  equal  to  the 
pore  size  a^,  that  is  to  say  Poiseuille  flow  breaks  down  when 


u  =2v/d^a 
c  f  p 


This  criterion  determines  if  the  viscositv  correction  factor  is  required. 


(b)  Permeability. 


The  permeabilty  k  of  a  porous  medium  is  related  to 


the  size  and  shape  of  the  pores  in  the  medium.  It  may  be  measured 
directlv,  calculated  from  measurements  of  the  flow  resistance,  or  esti¬ 
mated  with  one  of  several  commonly  used  empirical  equations.  When 

57  58 

experiments  are  carried  out  in  the  ocean,  the  Krumbein-Monk  or  Hazen 


relations  are  often  used;  these  equations  estimate  the  permeablity  of  a 
sediment  from  its  mean  grain  size.  Care  should  be  exercised  in  the  use 
of  these  equations  as  they  are  onlv  valid  in  regions  of  the  ocean  which 
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the  oscillatory  flow  of  the  viscous  fluid  through  the  pores  of  the 

medium.  The  effects  of  oscillatory  flow  were  first  studied  for  a  closed 

53 

channel  by  Kirchoff  in  1868,  oscillatory  viscosity  having  first  been 

54 

investigated  bv  Stokes  in  connection  with  the  oscillation  of  pendulums 
a  few  years  previously. 

For  porous  media,  Biot  assumed  that  at  low  frequency  the  flow  of  the 

12  55 

fluid  is  Poiseuille  flow  which  obeys  Darcy's  law.  ’  As  the  frequency 
of  oscillation  is  increased,  inertial  forces  in  the  flow  field  become 
comparable  with  the  viscous  forces  and  Darcy's  law  no  longer  adequately 
describes  the  flow;  eventually  the  flow  becomes  turbulent.  Biot  calcu¬ 
lated  a  complex  viscosity  correction  factor  whose  product  with  the  fluid 

13 

viscosity  v  accounts  for  the  deviation  from  Poiseuille  flow. 

The  complex  viscosity  correction  factor  was  defined  by  Biot  to  be: 

F(W)=( 1/4 )WT(W) /(l-2T(W)/iW)  (20) 

where , 

T(W)=[ber' (W)+i  bei ' (W) ] / [ ber(W)+i  bei(W)l  (21) 

and , 

W=a  (ud,/v) ' (22) 

P  f 

The  pore  size  parameter  a  is  usually  taken  to  be  one-sixth  to  one- 

P 

seventh  of  the  mean  size  of  the  individual  particles  when  the  Bint  model 
is  applied  to  sediments. 


»■; »  ‘  ■  1  1  17 .  i  n  ,k  "i 1  ■.  i  i  ■  v^r> 
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0.25.  The  values  of  A  and  n  for  a  wide  variety  of  materials,  may  be 
found  in  the  literature  (see  for  example  Ref.  50). 

(c)  Tortuosity  in  real  and  laboratory  sediments.  As  has  already  been 

stated,  tortuosity  has  been  the  most  widely  investigated  parameter  in  the 
Biot  model;  in  addition  to  its  relationships  with  formation  factor  and 
added  mass,  tortuosity  has  also  been  related  to  the  index  of  refraction 
of  fourth  sound  in  a  consolidated  (fused)  matrix  saturated  with  liquid 
Helium.  *  Experimentally,  the  tortuosities  measured  bv  the  electrical 
resistivity  method  and  tortuosities  measured  by  the  acoustic  index  of 

4 

refraction  of  superfluid  He  have  been  found  to  be  in  agreement  with  each 

. .  52 

other. 

For  real  and  artificial  sediments,  the  most  convenient  quantity  to 
measure  is  the  low-frequency  resistivity  formation  factor.  The 

A  Q 

appropriate  measurement  techniques  have  been  described  by  Wyllie  and 
45 

Sweet.  If  measurements  of  this  quantity  are  not  practical,  the 

41 

constants  required  for  use  of  the  Archie  equation  are  tabulated  for  a 

wide  variety  of  materials.50  For  closely  sized  spheres,  the  Slawinski 
49 

equation  can  be  used  to  compute  the  formation  factor  with  confidence. 

Use  of  the  Maxwell  expression  should  be  avoided  as  it  underestimates  the 

42-44 

formation  factor  for  spheres  in  contact. 

3.  Viscosity  and  permeability 

(a)  Viscosity.  If  the  fluid  independent  frame  moduli  frame  are  real 

quantities,  the  only  dissipation  mechanism  contained  in  the  Biot  model  is 


14 


F=(3-f)/2f. 


Measurements  performed  by  Boyce  in  the  Bering  sea  of  electrical 


conductivity  showed  an  electrical  conductivity  17%  below  that  predicted 

42 

by  the  Maxwell  expression;  this  experience  is  not  atypical.  Boyce 

derived  some  empirical  equations  which  fitted  his  data  better  than  the 

42  45 

Maxwell  expression,  as  have  many  other  investigator'".  *  The  expres- 

41 

sions  in  most  common  use  are  Archie's  Law  and  developments  of  it,  such 


as  the  Hutrule  formula  used  by  Schlumberger  in  well  loggine. 


It  is  not  surprising  that  the  resistivity  formation  factor  is  always 

underestimated  by  the  Maxwell  expression.  A  comparison  of  it  with 

46  47 

Berryman's  calculation  of  the  Hashin-Shtrikman  lower  bound  on  elec¬ 
trical  conductivity  shows  that  the  Maxwell  expression  is  identical  with 


this  lower  bound  on  F. 


Archie's  law  has  the  form, 


F=Af n, 


where  A  and  n  are  empirical  constants.  For  an  aggregate  of  unconsoli- 

48 

dated  spheres,  Wyllie  and  Gregory  found  that  the  formation  factor  could 

be  accurately  calculated  by  setting  A  =  1,  and  n  =  -1.3,  for  porosities 

49 

between  0. I  and  0.56.  They  further  found  that  the  Slawinski  equation, 


F=(l. 3219-0. 3219f)/f, 


fitted  all  available  data  for  sphere  packs  with  a  porosity  greater  than 


13 


different  values  of  the  coefficient,  but  added  mass  is  difficult  to 

calculate  for  a  sediment  matrix  because  it  is,  in  general,  a  tensor 
3  6“38 

quantity.  The  use  of  resistivity  formation  factors  seems  to  be  the 

most  promising  approach  for  ocean  acoustics  of  sediments. 

40 

(b)  Resistivity  formation  factor.  It  has  been  shown  by  Brown  that 

the  formation  factor  for  electrical  resistivity  and  tortuosity  are 
related  by 


X=Ff,  (15) 

where  the  resistivity  formation  factor  for  a  porous  medium  was  defined  by 
41 

Archie  as  the  ratio  of  the  resistivity  of  the  porous  sample  saturated 

with  a  conducting  fluid,  Rq,  to  the  resistivity  of  the  conducting  fluid 

in  bulk  R  ; 

w 

F=R  /R  .  (16) 

o  w 

This  equation  holds  only  if  the  matrix  material  is  an  insulator  and  the 
saturating  fluid  is  a  conductor.  The  resistivity  formation  factor  mea¬ 
sures  the  influence  of  the  pore  structure  on  the  conductivitv  of  the 
medium,  by  definition  it  is  always  greater  than  one. 

Little  work  has  been  done  on  the  electrical  conductivitv  of  ocean 

42  43 

floor  sediments.  It  was  long  expected  that  the  Maxwell  expression 

for  the  conductivity  of  an  assembly  of  spheres  embedded  in  a  conductor 

could  be  adapted.  The  Maxwell  expression  leads  to  the  following  eauation 

for  the  resistivity  formation  factor: 


t 


acceleration)  of  the  paddle  Is  greatly  increased  by  the  presence  of  the 
water.  This  increased  inertia  is  known  as  the  virtual  mass  of  the 
paddle;  the  difference  between  the  real  and  virtual  mass  of  the  paddle  is 
the  added,  or  induced  mass.  Thus,  the  added  mass  of  a  body  immersed  in  a 
fluid  is  the  difference  between  the  inertia  the  body  posseses  while  it  is 
in  the  fluid,  and  the  inertia  it  posseses  in  a  vacuum.  (The  preceeding 
illustration  is  due  to  Birkhoff.  ) 

It  is  customary  to  define  a  coefficient  of  added  mass  by  the  ratio, 

k  =  (added  mass  of  the  body)/(mass  of  fluid  displaced).  (13) 

The  added  mass  coefficient  is  a  function  only  of  the  shape  of  the  bodv; 

it  is  independent  of  the  size  of  the  body  and  its  acceleration,  and 

independent  of  the  density  of  the  fluid  in  which  the  body  is  immersed.  A 

sphere  oscillating  in  a  fluid  has  an  added  mass  coefficient  equal  to  one- 

38 

half.  Chapter  VI  of  Lamb's  treatise  on  hydrodynamics  is  devoted  to  the 

calculation  of  the  added  mass  of  nonspherical  bodies. 

39 

For  a  porous  medium,  Berryman  has  shown  that  the  added  mass  of  the 
frame  is  related  to  the  tortuosity  by 

X=l+k( f-1-l ) ,  (14) 

where  the  induced  mass  per  unit  volume  of  the  medium  is  eaual  to  kd^. 

For  an  aggregate  of  spheres,  the  added  mass  coefficient  for  an  isolated 

sphere  (1/2)  is  often  used.  This  yields  a  value  for  the  tortuosity  which 

is  too  high,  as  the  individual  spheres  are  not  as  free  to  move  as  a 

33 

single  sphere.  Some  authors  have  tried  to  compensate  for  this  by  usine 


compressibility  of  the  gas,  and  to  viscous  drag  at  the  pore  walls,  which, 
as  discussed  below,  is  the  only  attenuation  mechanism  originally  included 
in  the  Biot  model. 

2.  The  fluid-solid  coupling  factor  or  the  tortuosity 

The  density  parameter  which  appears  in  Biot's  coupled  differential 
equations 


d  -Xdf/f,  (6) 

where  X  is  the  fluid-solid  coupling  factor,  or  the  electrical  tortuosity 
of  the  frame,  is  the  most  investigated  -  and  best  understood  -  new 
parameter  introduced  into  the  acoustics  of  porous  media  by  the  Biot 
model.  X  is  often  referred  to  as  the  tortuosity  of  the  frame,  but  it 
should  not  be  confused  with  the  hydraulic  tortuosity  which  is  used  in  the 
Kozeny-Carman^  permeability  model.  The  hydraulic  and  electric  tortuosi¬ 
ties  are  related  to  each  other,  but  they  are  not  equal.  In  this 
dissertation,  when  tortuosity  is  referred  to,  electric  tortuosity  is 
always  implied.  Tortuosity  has  been  found  to  be  related  to  the  added 
mass  of  a  body  oscillating  in  a  fluid,  and  to  the  low-f requencv 
resistivity  formation  factor  of  the  matrix. 

(a)  Added  (or  induced)  mass.  Added  mass  has  been  discussed 

36  37 

extensively  by  Birkhoff  and  by  Landau  and  Lifshitz.  The  effects  of 
added  mass  are  often  experienced  in  the  macroscopic  world.  For  example, 


when  paddling  a  canoe, 


the  apparent  inertia  (i.e. 


resistance  to 


10 


K,  is  the  bulk  modulus  of  the  free  draining  frame,  and  G  is  the  shear 
b 

modulus  of  the  free  draining  frame.  In  specifying  a  free  draining  frame, 
we  mean  that  the  measurements  of  these  moduli  should  be  made  in  a  svstem 
which  allows  the  pore  fluid  to  escape  freely  as  the  sediment  is 


compressed. 


In  practice,  this  means  that  these  moduli  are  usually 


defined  in  terms  of  a  different  composite;  they  are  defined  for  a  matrix 

from  which  the  pore  fluid  has  been  removed.  The  moduli  K,  and  G  are 

b 

independent  of  the  pore  fluid  that  the  matrix  contains. 

The  dynamic  parts  of  the  moduli,  and  G,  are  then  defined  in  terms 

of  the  shear  and  compressional  wave  velocities  through  the  drained 

„  .  32 
matrix 


Vsh  =  {G/[(l-f)ds)} 


V  =  { [K  +(4/3)G]  /  [  (  1-f  )d  }}l,\ 
CO  s 


For  a  liquid  K,  and  G  are  in  general  complex  quantities.  The  real  parts 
b 

of  and  G  are  defined  through  Eqs.  (11)  and  (12);  their  imaginarv  Darts 

may  be  defined  through  measurements  of  the  logarithmic  decrement  of  a 

21  33™3A 

drained  frame,  or  from  empirical  relationships.  and  G  are 

allowed  to  be  complex  for  a  liquid  saturated  porous  material  in  order  to 

account  for  the  various  dissipation  mechanisms  which  occur  at  grain 

boundaries  (i.e.  friction,  relaxation  etc.).  If  the  pore  fluid  is  a  gas, 

then  the  frame  moduli  are  assumed  to  be  real  Quantities,  and  K.  and  G  are 

b 

completely  defined  by  Eqs.  (11)  and  (12);  for  a  gas  saturated  porous 
medium  all  the  attenuation  is  assumed  to  be  due  to  the  complex 


1.  The  Elastic  Constants 


H,  C,  and  M  are  elastic  constants,  which  are  in  general  complex.  The 

imaginary  parts  of  H,  C,  and  M  are  relevant  to  sound  attenuation  in  the 

medium.  H,  C,  and  M  may  be  defined  in  terms  of  more  familiar 
20  28-30 

moduli,  ’  some  real  and  some  complex: 

H“(Kr-Kb)2/(D-Kb)  +  +  4G/3  (7) 

C-K  (K  -K.  )/(D-K.  )  (8) 

r  r  b  b 

M=K  2/(D-K.  )  (9) 

r  D 

where , 

D-K  [l+f(K  /Kf-l)J  (10) 

The  most  lucid  description  of  the  "Gedanken"  experiments  through 

which  these  quantities  are  defined  is  given  by  Stoll  in  Ref.  20.  It  will 

be  sufficient  for  our  purposes  to  define  the  terms  appearing  in  Eqs.  (7)- 

(10),  and  to  explore  their  physical  content. 

is  the  bulk  modulus  of  the  individual  particles  which  make  up  the 

sediment  matrix,  this  quantity  is  always  a  real  quantity.  It  is  assumed 

that  this  bulk  modulus  is  for  isothermal  conditions.  K^.  is  the  bulk 

modulus  of  the  pore  fluid.  The  bulk  modulus  of  the  pore  fluid  is  assumed 

to  be  a  real  quantity  if  the  pore  fluid  is  a  liquid;  if  the  pore  fluid 

3 1 

is  a  gas,  its  bulk  modulus  is  sometimes  regarded  as  a  complex  quantity 
to  account  for  sound  attenuation. 


> 

• 

8 

•  '  -  _*  •  * 

i 

A  change  of  variables  is  normally  made  when  the  Biot  model  is  applied 

• 

to  marine  sediments,  so  that  the  equations  appear  in  terms 

of  the  dilita- 

tion  of  the  frame  e,  and  the  relative  dilitation  between 

the  frame 

and 

i 

the  fluid  r,  rather  than  in  terms  of  the  dilitation  in  the  fluid  E 

and 

• 

the  dilitation  in  the  solid  e, 

r 

r=f  div(u-U)=f (e-E). 

(3) 

i 

• 

It  is  this  set  of  transformed  equations  which  will  be  considered  in 

the 

following: 

» 

The  coupled  equations  describing  the  propagation  of 

compressional 

• 

waves  are: 

V2(HE-Cr)=D2(de~dfr), 

/AN 

1 

) 

- 

•  , 

and 

*  *  ^4  * 

K 

V2(Ce-Mr)=D2(dce-d  r)  -  (v/k)F(W)D  r, 

t  L  C  t 

(5) 

•  , 

where 

» 

d  -Xd-/f. 
c  f 

(6) 

- 

• 

-  4 

denotes  the  kth  partial  derivative  with  respect  to  the 

variable  i 

• 

The  differential  equation  for  shear  waves  will  not  be 

presented 

here , 

» 

• 

but  the  dispersion  equation  for  shear  waves  is  given  in  section  II 

.  B.  4. 

The  various  terms  appearing  in  these  equations  will  be 

defined  , 

and 

discussed  in  some  detail  below. 

» 

- 

• 

V .•  ,‘.’l 

» 

• 
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to  be  homogeneous  and  isotropic.  The  size  of  the  pores  is  required  to  be 
small  compared  to  the  wavelength  of  any  disturbance.  If  this  condition 
is  not  met  scattering  occurs  and  the  model  is  not  applicable. 

The  equations  of  motion  for  a  poro-elastic  medium  are  derived  from 
the  stress-strain  equations.  By  requiring  the  theory  to  be  linear,  the 
potential  energy  can  be  expressed  as  a  quadratic  function  of  the  strain 
components,  and  the  kinetic  energy  can  be  expressed  as  a  quadratic 
function  of  the  fluid  and  solid  velocities.  The  average  stress  acting  on 
an  arbitrary  volume  element  is  defined  to  be  the  sum  of  the  forces  acting 
on  the  fluid  and  solid  components  divided  by  the  area  of  the  volume 
element.  Strain  is  defined  in  terms  of  the  displacements  of  the  fluid  U_, 
and  the  solid  u.  A  Lagrangian  analysis  is  then  performed  resulting  in  a 
pair  of  simultaneous  coupled  differential  equations  in  terms  of  the  fluid 
and  solid  displacements.  These  equations  are  then  separated  into  a  nair 
of  simultaneous  coupled  differential  equations  containing  only 
dilitations,  this  pair  describes  the  propagation  of  compressional  waves; 
and  a  pair  of  simultaneous  coupled  differential  equations  containing  onlv 
rotations,  describing  the  propagation  of  shear  waves. 

In  these  equations,  the  aggregate  density,  d,  of  the  fluid-saturated 
porous  medium  is  related  to  the  density  of  the  fluid,  d^,  and  the  density 
of  the  grains,  dg,  bv  the  equation, 

d-(l-f)d  +  fdf,  (2) 
s  t 

where  f  is  the  porosity  of  the  matrix. 


6 


applicable  as  the  Biot  model;  the  Biot  model  will  reduce  to  the  other 

models  under  the  restrictive  conditions  to  which  they  apply.  Biot  has, 

in  subsequent  years,  further  generalized  the  model,  he  has  also  simpli- 

1 5“  1 8 

fied  the  derivations  of  the  fundamental  equations.  The  model  was 

19 

not  applied  to  sediments  until  1969,  when  it  was  done  by  Stoll  ;  the 

version  of  the  theory  presented  in  this  dissertation  is  normally  known  as 

the  Biot-Stoll  model,  although  most  of  the  necessary  developments  may  in 

28 

fact  be  found  in  the  earlier  work  of  Biot. 

The  principal  prediction  of  the  Biot  model  is  that  a  fluid  saturated 
porous  medium  will  support  three  body  waves.  Two  compressional  waves, 
and  one  shear  wave  are  predicted.  The  shear  wave  is  essentially  similar 
to  shear  waves  which  propagate  in  linear  elastic  solids;  however  it  is 
modified  by  the  presence  of  the  fluid.  The  compressional  waves  are 
usually  known  as  the  fast  wave  (corresponding  to  in-phase  motion  of  the 
fluid  and  the  frame),  and  the  slow  wave  (corresponding  to  out-of-phase 
motion  of  the  fluid  and  the  frame).  These  waves  are  also  known  as  waves 
of  types  I  and  II,  respectively.  The  fast  wave  is  again  essentiallv  the 
compressional  wave  of  a  linear  elastic  solid,  the  slow  wave  is  a  new  wave 
which  is  charactersitic  of  wave  propagation  in  fluid  saturated  porous 
media.  The  properties  of  these  waves  will  be  further  discussed  in  sec¬ 
tion  II. B. 4,  where  the  dispersion  equations  of  three  bodv  waves  are 
presented. 

Biot  treated  a  fluid-saturated  porous  medium  as  a  solid  matri..  filled 
with  arbitrarily  shaped  pores.  The  distribution  of  the  pores  is  assumed 


5 


the  frequency  of  the  sound  wave  measured  i  kilohertz,  and  n  is  the 

exponent  of  frequency  taken  to  be  equal  to  one  in  Hamilton's  model. 

The  attenuation  of  sound  in  a  particular  sediment  may  then  be 

predicted,  if  the  sediment  type  is  known,  that  is  to  say  that  the 

density,  porosity,  or  particle  size  of  the  sediment  is  known.  Values  of 

the  proportionality  constant,  K,  are  tabulated  and  presented  graphically 

as  functions  of  the  density,  porosity,  or  mean  particle  size  in  many  of 
6_8 

Hamilton's  papers  for  a  wide  variety  of  sediments. 

The  Hamilton  model  assumes  that  there  is  no  significant  velocity 
dispersion  for  compressional  waves  in  sediments  from  frequencies  of  a  few 
hertz  to  frequencies  of  several  megahertz.  Recent  data  casts  doubt  on 
the  validity  of  this  assumption,1^  11  and  as  this  assumption  is  the  basis 
of  the  Hamilton  model,  another  model  is  clearly  required. 

B.  The  Biot  Model 

The  Biot  model,  which  is  now  the  most  widely  accented  model  for  sound 

propagation  in  fluid-saturated  porous  media,  was  originallv  proposed  bv 

12-13 

Biot  in  1956,  in  two  papers.  The  first  paper  dealt  with  wave  propa¬ 

gation  at  low  frequencies,  and  the  second  dealt  with  wave  propagation  at 
high  frequencies.  The  terms  low  and  high  frequency  as  applied  to  the  Biot 
model  will  be  defined  below  in  section  II.B.3.a.  The  development  of  the 

model  used  some  results  from,  and  the  notation  of,  Biot's  earlier  work  on 

22 

soil  consolidation  problems. 

Similar  theories  had  been  developed  in  earlier  years  bv  other 
23-27 


authors 


but  none  of  these  other  theories  was  as  generallv 


II.  THEORETICAL  MODELS 


In  this  chapter  two  widely  used  models  of  sound  propagation  in 
sediments  are  reviewed.  The  models  are  known  as  the  Hamilton  and  Biot 
models,  although  the  latter  is  sometimes  called  the  Biot-Stoll  model. 
The  Biot  model  is  presented  in  somewhat  greater  detail,  both  because  of 
its  greater  complexity,  and  because  the  development  of  this  model  took 
place  in  a  long  series  of  papers,  and  no  adequate  review  of  the  model 
exists.  Although  the  primary  concern  of  this  dissertation  is  the  atten¬ 
uation  of  compressional  (sound)  waves,  some  attention  will  be  given  to 
the  attenuation  of  shear  waves. 

Throughout  this  dissertation  attenuations  denoted  by  the  svmbol  a  are 
measured  in  dB/m,  while  those  denoted  by  the  symbol  a  are  measured  in 
nepers/m. 


A.  The  Hamilton  Model 

Hamilton  has  been  investigating  the  acoustic  properties  of  marine 
sediments  since  the  early  1950s. ^  ^  The  Hamilton  model  predicts  that  the 
attenuation  of  compressional  waves  in  sediments  is  proportional  to  the 
first  power  of  frequency.  That  is  in  the  equation, 

a=Kfn,  (1) 

a,  is  the  attenuation  of  compressional  waves  measured  in  dB/m,  K  is  a 
proportionality  constant  dependent  on  the  sediment  classification,  f  is 


have  sediment  types  similar  to  the  area  in  which  the  measurements  were 

made  from  which  the  relations  were  extracted.  They  are  not  general  and 

if  measurements  are  possible  they  arc  to  be  preferred.  In  laboratory 

experiments  using  closely-sized  spheres  the  Kozenv-Carman  relation  mav  be 

35 

used  with  confidence. 

4.  Solution  of  Biot's  equations 

For  an  infinite  medium,  the  dispersion  equation  for  compressional 

waves  in  a  porous  medium  may  be  readily  found  by  assuming  that  Biot's 

equations  admit  plane  wave  solutions  and  substituting  these  solutions 

into  the  coupled  differential  equations.  For  wave  propagation  in  bounded 

media,  appropriate  boundary  conditions  have  been  derived  by  Deresiewicz 
59 

and  Skalak.  Deresiewicz  has  investigated  the  properties  of  several 

different  types  of  boundaries  in  a  long  series  of  papers. ^  ^ 

Biot's  equations  have  been  solved  in  cylindrical  seometrv  bv 
Gardner^  and  Berryman. ^  Gardner's  solution  is  valid  only  for  very  low 
frequencies  while  Berryman’s  solution  is  valid  for  the  same  range  of 
frequencies  for  which  the  Biot  model  is  valid.  Both  of  these  investiga¬ 
tions  dealt  with  the  propagation  of  extensional  waves  in  consolidated 
porous  media;  it  is  not  immediately  obvious  that  they  are  directlv 
relevant  to  the  experiments  performed  in  this  investigation  as  plane 
waves  are  generated.  It  is  expected  that  if  the  cylindrical  boundaries 
do  have  any  effect  on  the  attenuation  that  it  could,  in  the  first 
approximation,  be  accounted  for  by  the  classical  solution  of  Pochhammer 


and  Chree 


73-76 


for  wave  propagation  in  elastic  rods 


19 


To  obtain  the  dispersion  equation  for  congressional  waves,  we  assume 
plane  wave  solutions  of  the  form: 


e=A^exp  [  i(u)t-lx)  ] 


r=A2exp[i(dJt-lx)  ] . 


Substitution  of  these  solutions  into  Biot's  equations  leads  to  the  deter¬ 


minants!  frequency  condition, 


2  2 
HI  -da 

2  2 
Cl  -dfw 


2  2 
dfw  -Cl 

d  u)2-Ml2-iioF(U)v/k 
c 


The  roots  of  this  equation  give  the  attenuation  1  ,  and  the  phase 

1m 

velocity  w/lRe  as  functions  of  frequency  for  compressional  waves  of  the 
first  and  second  kinds.  A  similar  procedure  is  used  to  obtain  the 
dispersion  equation  for  shear  waves. 

The  determinant  is  usually  solved  numerically,  however  the  complex 
wavenumbers  for  shear  and  compressional  waves  will  be  given  here  for 
completeness 

u2{(Md-2Cd  +Hq)-((Md-2Cd  +Hq)-4(C2-HM)(d  2-dq)] l/2} 

12=  - - - ^ - - -  (28) 

r  2(HM-C  ) 


2  221/2 
,  ui  {(Md-2Cdf+Hq)+[(Md-2Cd,+Hq)-4(C  -HM)(d  -da) )  }, 

if-  - £ - - - - -  (29) 

2(HM-C") 


For  shear  waves, 


1sh2='1)2(d~df2/q)/G’  O0) 

where , 

q=d  -iF(W) v/k  .  (31) 

c 

In  these  equations  f,s,sh  on  left  hand  sides  denotes  fast  compressional, 

slow  compressional ,  and  shear  waves  respectively. 

All  three  of  the  body  waves  are  attenuated  and  dispersive.  At  low 

frequencies  the  slow  compressional  wave  is  not  a  propagational  wave.  The 

relative  effects  of  the  various  parameters  contained  in  these  equations 

1 2-21 

on  wave  velocity  and  attenuation  have  been  widely  reported  and  will 

not  be  repeated  here. 

5.  Conclusions 

It  should  not  be  too  surprising  that  there  are  three  propagational 
body  waves  in  porous  media  as  Biot  treated  the  fluid  and  solid 
constituents  of  the  medium  on  an  equal  footing.  A  fluid  will  support  a 
compressional  wave,  and  the  matrix  material  will  support  a  compressional 
wave  and  a  shear  wave.  All  three  waves  have  been  detected 
experimentally,  although  the  compressional  wave  of  the  second  kind  was 
elusive,  and  was  only  detected  in  recent  years  by  Plona.^ 

Biot's  equations  have  been  successfully  applied  to  a  varietv  of 
physical  systems,  including  the  propagation  of  fourth  sound  in  a 


superfluid/superleak  system  ,  pressure  diffusion  through  porous 

81 "8 2  83 

media,  ,  the  elastodynamics  of  gels  ,  and  slow  waves  and  the 

consolidation  transition. ^ 

The  main  difficulty  in  applying  the  model  to  sediments  is  in  the 

measurement  or  calculation  of  the  fluid  independent  frame  moduli  K,  and 

b 
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G.  For  a  consolidated  matrix,  Berryman  has  developed  self-consis¬ 

tent  schemes  for  their  calculation,  but  these  are  unsuitable  for  the 
unconsolidated  or  partially  consolidated  matrix  typical  of  surficial 
ocean  sediments.  The  most  commonly  used  method  is  to  assume  a  value  for 
Poisson's  ratio  for  the  sediment  and  to  vary  and  G  until  the  Biot 
model  produces  the  measured  phase  velocities.  Future  theoretical  work 
will  have  to  concentrate  on,  either  the  calculation  of  H,C,  and  M 
directly,  or  a  satisfactory  method  of  calculating  the  frame  moduli  for 


unconsolidated  media 


III.  The  Sediment 


A.  Real  and  Laboratory  Sediments 


This  investigation  reports  measurements  performed  on  laboratory 
(artificial)  sediments.  The  selection  of  an  artificial  sediment  was  not 
an  easy  one,  as  one  of  the  long-range  aims  of  this  investigation  is  to 
provide  a  better  understanding  of  the  relationship  between  theoretical 
models  of  attenuation  and  in-situ  attenuation  measurements.  However,  the 
obvious  short-term  advantages  of  using  a  real  sediment  were  outweighed  by 
the  requirement  of  a  homogeneous,  isotropic  sample.  The  Biot  model 
describes  particles  and  pore  spaces  which  are  uniformly  sized,  although 
it  has  been  applied  to  sediments  containing  a  wide  range  of  sizes.  For 
this  reason,  a  well  characterized  sample  of  closely  sized  spherical  glass 
beads  was  chosen  for  the  initial  measurements.  Real  and  laboratory 


sediments  differ  in  some  important  ways,  we  shall  now  consider  some  of 
these  differences. 

$7  p 

In  an  in-situ  sediment,  the  density,  porosity,  and  bulk  and 

89 

shear  moduli  of  the  sediment  are  all  functions  of  deoth  within  the 
sediment,  and  of  the  overburden  pressure.  The  velocity  and  attenuation 


of  shear  and  compressional  waves  are  also  functions  of  depth.  There  is  a 

7  38—89 

lot  of  data  available  on  velocity  gradients  in  sediments,  ’  but 


little  available  on  the  variation  of  attenuation  with  depth. 


10,88 


If  the 


porosity  of  a  sediment  is  not  a  function  of  depth,  the  sediment  is 

90 

described  as  a  non-compactive  sediment.  Most  in-situ  sediments  are 


23 


compactive,  but  an  artificial  sediment  may  be  constructed  so  as  to  be 
non-compacti ve . 

In  this  investigation  it  would  clearly  be  advantageous  to  construct  a 
non-compactive  sediment.  In  a  non-compacti ve  sediment  the  densitv  and 
frame  bulk  modulus  of  the  sediment  are  not  functions  of  depth;  the  rela¬ 
tionship  between  aggregate  density  and  porosity  was  given  in  the  last 

chapter;  the  bulk  modulus  of  the  frame  is  also  a  function  of  porosity  in 

90 

the  first  approximation.  The  shear  modulus  of  the  frame  is  still  a 


function 

of  depth  in  a  non-compactive 

sediment;  this  is 

because 

the 

deeper  a 

particle  is  within  the  body 

of  the  sediment, 

the 

more 

its 

movement 

is  restricted  by  the  weight 

of  particles  above 

it. 

90 

Thus , 

given  the 

functional  dependence  of  the 

complex  wave  numbers 

of 

shear 

and 

compressional  waves  on  the  shear  modulus,  the  velocity  and  attenuation  of 
shear  and  compressional  waves  should  be  expected  to  remain  functions  of 
depth.  This  effect  is  expected  to  be  small  in  the  case  of  comoressional 
waves  as  the  shear  modulus  is  usually  at  least  an  order  of  magnitude 
smaller  than  the  other  material  moduli  in  an  unconsolidated  sediment. 


There  are 

some 

other  differences  between 

real  and  laboratory 

sedi- 

ments  which 

should 

be  noted. 

In  natural 

sedimentary  materials 

the 

sphericity  of 

the 

particles 

decreases  with 

decreasing  size  of 

the 

particles,  so  small  spherical  glass  beads  are  not  a  good  representation 
of  similarly  sized  particles  in  a  real  sediment.  In  finely  sized  sedi¬ 
ments  the  particles  are  often  angular.  Spherical  glass  beads  will  tend 
to  have  a  smaller  shear  modulus  as  their  smooth  surfaces  allow  them  to 


move  more  easily.  Compressional  wave  velocities  are  higher  in  laboratory 
sediments  than  they  are  in  real  sediments  because  the  bulk  modulus  of 
glass  is  higher  than  the  bulk  modulus  of  quartz. 

B.  Physical  Properties  of  the  Sediment 

1.  The  glass  beads 

The  glass  beads  used  in  this  experiment  were  obtained  from  the  Ferro 

Corporation  in  Jackson,  Mississippi.  The  manufacturer  describes  the 

beads  as  Class  III  Glas-Shot  (sic),  they  are  sold  for  use  as  abrasives. 

The  sales  literature  gives  the  following  description  of  the  beads: 

Glas-Shot  (beads)  are  spheroidal  in  shape  containing  not  more 
than  fifteen  (15)  percent  irregularly  shaped  particles.  They 
are  reasonably  free  of  sharp  angular  particles,  the  quantity 
of  which  does  not  exceed  three  (3)  percent,  oarticles  showing 
milkiness  or  surface  scoring  or  scratching,  the  quantity  of 
which  does  not  exceed  two  (2)  percent,  and  foreign  matter,  the 
quantity  of  which  does  not  exceed  one  (1)  percent  including  iron 
particles,  the  quantity  of  which  does  not  exceed  one-half  (1/2) 
percent. 

The  beads  selected  are  designated  as  type  MS-L  by  the  manufacturer.  The 
beads  range  in  size  between  44  and  88  microns.  A  summary  of  the 
manufacturers  specifications  of  the  bulk  properties  of  the  glass  is  given 
in  Table  1.  The  glass  is  a  soda-lime-silica  glass.  A  glass  of  this 
composition  is  usually  called  a  soda-lime  glass  and  will  be  so  called 
throughout  this  dissertation. 

There  is  some  ambiguity  over  which  modulus  of  elasticity  is  provided 

by  the  manufacturer,  as  only  modulus  of  elasticity  is  specified.  Consul- 
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tation  with  a  manual  of  glass  properties  suggests  that  the  modulus  of 
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TABLE  1 


Manufacturer's  Specifications  of  Glass  Beads 


Sound  Velocity  (glass  in  bulk) 
Modulus  of  Elasticity 
Density 
Thermal  Conductivity 
Coefficient  of  Thermal  Expansion 
Mean  Specific  Heat 


5800  ms 

9.66  x  10^  psi 
2.5  g  cm 

2.5  x  10  3  Cal  s  *  cm  3 
8.47  x  10~6 

0.27  Cal/g/°C  (40-800  °C) 


Properties  of  Glas-Shot  Beads 


Bulk  Specific  gravity  1.5  (average) 

1.55  (coarse  sizes) 
1.45  (fine  sizes) 

Diameter  of  MS-L  Beads  44-88  microns 
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elasticity  quoted  by  the  manufacturer  is  a  typical  value  of  Young's 
modulus  for  a  soda-lime  glass. 


2.  Porosity  and  laboratory  sediments 


The  density  and  porosity  of  a  regular  packing  of  identical  spheres 

does  not  change  with  grain  diameter,  but  a  regular  packing  of  identical 

spheres  does  not  have  regularly  sized  pores,  and  its  acoustic  properties 

are  aeolotropic.  The  porosity  of  a  randor,  packing  of  uniformly  sized 

spheres  is  0.38.  Variations  from  this  ideal  porosity  are  caused  by 
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variations  in  grain  size  and  shape,  or  in  the  packing. 
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Dullien  defines  four  modes  of  random  packing.  We  need  not  concern 
ourselves  with  the  details  of  how  these  oackings  may  be  achieved,  but 
should  note  that  variations  in  porosity  from  0.36  to  0.44  are  possible 
for  random  packings  of  spheres. 

Some  measurements  of  porositv  were  carried  out  on  small  samples  in 
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the  laboratory  using  the  liquid  imbibtion  techniaue.  When  water  was 
added  to  the  beads  a  porosity  of  0.38  was  measured.  When  beads  were 
added  to  water  a  porosity  of  0.41  was  measured;  when  the  samples  were 
allowed  to  sit  overnight  the  measured  porosity  decreased  to  0.39.  The 
standard  deviation  in  all  these  porositv  measurements  was  0.01.  A 
porosity  of  0.40  is  used  in  all  subsequent  calculations. 

In  the  sound  tube,  beads  were  added  to  water.  Although  some  settling 
takes  place,  fewer  gas  bubbles  are  trapped  in  this  wav.  In  the  porosity 
measurements,  when  water  was  added  to  the  beads,  it  took  a  long  time  to 


permeate  the  column  and  the  sample  had  to  be  stirred  to  remove  some  of 


the  bubbles. 

The  container  walls  will  have  some  effect  on  the  packing  of  the 

beads,  and  therefore  the  porosity.  The  base  plane  of  the  Dack  will 

introduce  local  order  into  a  random  packing,  but  even  if  this  first  layer 

is  completely  regular,  the  second,  and  all  subsequent  lavers  will  be 

built  from  multiple  focal  points  and  the  pack  will  quickly  assume  the 

characteristics  of  a  random  pack.  Near  any  external  surface  there  will 

be  a  region  of  high  porosity  because  of  the  difference  in  the  radii  of 

curvature  of  the  beads  and  the  container  walls.  The  wall  effects  on  bulk 

porosity  are  negligible  if  the  diameter  of  the  container  is  more  than  ten 
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times  the  mean  diameter  of  the  particles  in  the  pack.  That  condition 
was  met  in  the  experiment  described  here. 

3.  Permeability 

Permeabilty  is  the  term  used  for  the  conductivity  of  a  porous  medium 

with  respect  to  permeation  by  a  Newtonian  fluid.  Permeability  was  not 

measured  directly  in  this  investigation;  instead  the  flow  resistance  of 
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the  medium  was  measured  using  a  Leonard  flow  resistance  apparatus 
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modified  as  suggested  by  Rudnick.  The  permeability  of  a  porous  medium 
is  equal  to  the  viscosity  of  the  permeating  fluid  divided  bv  the  flow 
resistance  of  the  medium.  This  method  measures  the  gas  permeabilitv 
rather  than  the  absolute  permeability.  The  absolute  permeabiltv  mav  be 
calculated  by  plotting  gas  permeability  versus  the  inverse  of  the  average 
pressure  across  the  porous  medium.  In  this  investigation,  onlv  the 


differential  pressure  across  the  sample  was  measured. 


The  measured  gas 


permeabilities  are  presented  in  table  2. 

The  permeability  of  a  porous  medium  consisting  of  uniformly  sized 

35 

spherical  particles  may  be  calculated  from  the  Kozeny-Carman  equation 

k=f3/[k  S  2(l-f2)]  (32) 

o  o 

where  S  is  the  specific  surface  of  the  particles,  and  k  is  an  empirical 
o  o 

constant.  For  uniform  spheres, 

SQ=6/d,  (33) 

where  d  is  the  diameter  of  the  particles.  The  constant  k  mav  be  is 

o 

equal  to  five  for  a  wide  variety  of  porous  media;  thus  the  permeabiltv  of 
the  sphere  pack  may  be  calculated  from: 

k=d2f3/[180(l-f2)].  (34) 

The  particles  used  in  this  investigation  have  diameters  between  44  and  88 

microns.  These  beads  have  been  used  previously  bv  Bell.^  He  found  that 

their  mean  diameter  was  62.5  microns.  Using  this  diameter,  the  Kozenv- 

-12  2 

Carman  equation  predicts  a  permeabiltv  of  3.86  x  10  m  at  a  porositv  of 
0.40;  this  compares  favorably  with  the  measurements  reported  in  table  2 
and  is  used  in  all  subsequent  calculations.  A  summary  of  the  physical 
pronerties  of  the  sediment  is  provided  in  table  3. 
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TABLE  2 


Permeability  and  Flow  Resistance  measured  in  Leonard's  Apparatus 


.f ferential 
Pressure 

Flow 

Resistance 

Permeability 

(Nm“2) 

(cgs  Rayls) 

12  2 
(10iZ  tn  ) 

- 

96.18 

5045 

3.58 

192.35 

4742 

3.82 

384. 7 1 

4378 

4.13 

480.88 

4242 

4.27 

673.24 

4440 

4.08 

961.76 

4348 

4.16 

— .1*11.111^..  ■>*.. .Jl  wfti'afT  ill  ^ 
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TABLE  3 


Summary  of  Sediment  Physical  Properties 


Density  of  Grains 

2500  kg/m3 

Density  of  Fluid 

998  kg/m3 

Bulk  Modulus  of  Grains 

3.83  x  1010  1 

Bulk  Modulus  of  Fluid 

2.0  x  109  Pa 

Viscosity  of  Fluid 

1.0  x  10”3  N 

Mean  Particle  Size 

62.5  x  10"6  ■ 

Permeability  of  Frame 

3.86  x  10”12 

Porositv  of  Frame 


40% 


IV.  EXPERIMENT 


The  velocity  of  compressional  sound  waves  in  sediments  is  more  than 

4-10 

1500  meters  ner  second;  this  means  that  laboratory  measurements  of 

the  low-frequency  attenuation  of  compressional  sound  waves  in  sediments 
require  a  large-scale  apparatus  so  that  the  length  of  the  sediment  sample 
is  longer  than  one  wavelength  of  the  compressional  wave  at  the  lowest 
desired  frequency. 

This  chapter  contains  a  description  of  the  experimental  apparatus 
used  to  perform  low-frequency  attenuation  measurements  and  details  of  the 
experimental  procedure  employed. 

A.  Experimental  Apparatus 


The  attenuation  measurements  reported  in  this  dissertation  were 
performed  in  a  vertically  mounted  sound  tube.  The  tube  is  made  of  alumi¬ 
num,  it  is  4.35  m  long,  and  has  an  inside  diameter  of  0.254  m  with  a  wall 

thickness  of  0.0095  m.  The  working  length  of  the  sound  tube  is  3.75  m, 

3 

which  enables  the  tube  to  hold  a  maximum  of  0.19  m  (50  gallons)  of 
sediment . 

The  sound  tube  is  mounted  on  three  wooden  supports,  and  is  surrounded 


by 

a  scaffold  to 

allow  easv  access  to  all 

parts 

of  the  tube. 

There 

are 

sui 

table 

svstems 

which  allow  water  to  be 

added 

to  or  removed 

f  rom 

Che 

svs 

tem. 

The  system  is  diagrammed  in  Fig. 

1. 
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It  is  permissible  to  treat  a  sediment  as  a  fluid  for  the  purpose  of 

letermining  wall  losses  in  a  clyindrical  tube,  provided  that  the  sediment 
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ihear  modulus  is  small  in  comparison  with  the  frame  bulk  modulus; 
:his  condition  is  met  in  an  unconsolidated  sediment  which  is  bv 
lefinition  a  sediment  of  low  rigidity.  Kirchhoff's  wall  loss  expression 
fas  used  to  correct  the  attenuation  data  presented  in  Fig.  5  for  the 
ittenuation  due  to  the  walls  of  the  sound  tube. 

In  a  sound  tube  of  radius  R  for  a  sound  wave  of  frequency  v 
propagating  in  a  medium  of  viscosity  v,  thermal  conductivity  K,  with  a 
ratio  of  specific  heats  Y,  and  an  average  fluid  pressure  p,  Kirchhoff 
.fall  losses  may  be  calculated  from: 


a  „  IfJL  +  Slzll  (JL))(2Vv ji/2 
w  R  yh  +  y  (vcv,i(  P  } 


where  a  is  the  sound  attenuation  due  to  the  walls  of  the  tube. 


specific  heat,  C,  of  a  sediment  may  be  calculated  from: 


(35) 


The 


dC  =  (l-f)d  C  +  fdcCc  (36) 

s  s  t  f 

where  f  is  the  frame  norosity,  and  d,  d^  and  d  are  the  densities  of  the 
aggregate  the  solid  and  the  fluid.  The  ratio  of  specific  heats  for  water 

3 

is  1.004.  The  specific  heat  of  water  (C^)  is  4184  J/m  °C  and  the 

3 

specific  heat  of  a  sode-lime  glass  (C  )  is  840  J/m  'C;  thus  using  the 

s 

expression  for  the  specific  heat  of  a  sediment  the  ratio  of  specific 

heats  for  the  sediment  is  calculated  to  be  1.002. 

The  thermal  conductivity  K  of  a  sediment  mav  be  estimated  from  its 
90 


water  content  using: 


1  2  3  456789  10 


Frequency  (kHz) 

Figure  5 

Attenuation  Versus  Frequency  Uncorrected  for  Tube  Absorption 
(Bars  indicate  scatter  between  two  measurements) 


V.  RESULTS  AND  DISCUSSION 


The  attenuation  of  compressional  waves  in  a  laboratory  sediment  was 
determined  from  the  received  signal  level  over  different  oath  lengths 
through  the  sediment.  Path  lengths  of  2.114  m,  2.206  m,  2.679  m,  and 
3.073  m  were  used.  The  measurements  at  2.114  m  and  2.206  m  were  maae 
with  a  different  transmitter  to  the  measurements  made  at  2.679  m  and 
3.073  m.  The  attenuation  measurements  reported  in  this  dissertation  were 
calculated  by  matching  the  gradients  of  received  signal  level  versus  oath 
length  for  the  two  transmitters  to  determine  the  average  attenuation.  No 
systematic  variation  in  the  gradients  produced  by  the  different  transmit¬ 
ters  was  noted.  The  resulting  attenuation  measurements  are  dotted 
versus  frequency  in  Fig.  5. 

When  sound  waves  propagate  in  a  cylindrical  tube  there  is  some 

attenuation  at  the  walls  produced  by  viscosity  and  thermal  conduction. 

The  theory  of  wall  losses  in  a  cylindrical  tube  containing  a  gas  was 
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worked  out  by  Kirchhoff  in  1968.  It  is  possible  to  demonstrate  that 
Kirchhoff's  expression  for  wall  losses  is  valid  for  a  cylinder  containing 
any  fluid. When  low-pressure  gases  are  used  it  is  necessarv  to  intro¬ 
duce  velocity  slip  and  temperature  jump  at  the  walls  into  Kirchhoff's 
solution; these  additional  terms  are  unnecessary  if  the  fluid  is  a 
water  saturated  porous  medium  because  the  mean  free  path  in  such  a  medium 
is  on  the  order  of  atomic  dimensions. 
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More  sediment  was  then  added  to  the  tube  using  the  methods  previously 
outlined  in  this  section. 

The  experimental  technique  used  frequently  caused  damage  to  the 
transmitter.  The  damage  occurred  for  two  reasons.  Glass  beads  stick  to 
the  walls  of  the  sound  tube.  These  glass  beads  are  highly  abrasive  and 
can  cause  significant  damage  to  the  transmitter.  The  seal  formed  bv  the 
Teflon  gasket  surrounding  the  transmitter  when  it  is  in  contact  with  the 
water  can  also  cause  problems.  On  one  occasion  it  caused  air  to  be 
forced  behind  the  Mylar  film  as  the  transmitter  was  removed  from  the  the 
tube.  The  attenuation  measurements  reported  in  Chapter  V  were  made  with 
two  different  transducers;  none  of  the  transducers  survived  more  than 
three  immersions. 

Since  this  transmitter  was  designed  and  constructed,  better  designs 

for  underwater  solid  dielectric  transducers  have  appeared  in  the  litera- 
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ture.  Experience  gained  during  the  course  of  this  investigation 

suggests  that  it  might  be  advantageous  to  fill  the  tube  with  water  above 
the  level  of  the  sediment  and  to  correct  the  resulting  measurements  for 
the  (small)  attenuation  due  to  the  water. 


sealed  systems  was  then  evacuated.  The  pressure  in  the  system  was  main¬ 
tained  at  the  vapor  pressure  of  water  for  several  hours  to  allow  anv 
entrapped  gas  to  escape  from  the  surface  of  the  sediment.  The  vacuum  was 
then  released  and  the  length  of  the  sediment  sample  was  determined.  The 
level  of  water  above  the  glass  beads  was  usually  about  six  inches. 

After  the  sample  had  been  prepared,  the  transmitter  was  introduced 
into  the  sound  tube.  The  transmitter  was  allowed  to  sink  to  the  level  of 
the  water  under  its  own  weight.  A  continuous  signal  was  broadcast  by  the 
transmitter  so  that  by  monitoring  the  signal  received  by  the  hydrophone 
we  could  deterine  when  the  transmitter  was  in  contact  with  the  water.  No 
clear  signals  were  observed  until  the  transmitter  was  in  contact  with  the 
surface  of  the  water.  With  experience  it  became  simple  to  recognize  when 
the  transmitter  was  in  contact. 

When  the  transmitter  was  in  contact  with  the  water,  the  electronics 
were  adjusted  so  that  tone-bursts  were  transmitted  through  the  sediment 
and  the  received  signal  level  was  adequate.  The  same  amplifier  settings 
were  then  used  in  all  subsequent  experiment.  Compressional  waves  with 
frequencies  between  2  kHz  and  10  kHz  were  transmitted  through  the  sedi¬ 
ment.  The  peak-to-peak  output  voltage  of  the  received  signal  as 
amplified  by  the  conditioning  amplifier  and  the  PAR  pre-amnli f ier  and 
displayed  on  the  oscilloscope  was  recorded. 

In  order  to  determine  sound  attenuation,  the  oath  length  through  the 
sediment  must  be  varied.  After  the  completion  of  a  data  run  the 
transmitter  was  removed  from  the  sound  tube  using  a  block  and  tackle. 


4.  The  water  supply 


Water  is  added  to  or  removed  from  the  system  bv  a  pair  of  water 
valves  mounted  to  the  bottom  end-plate  of  the  sound  tube.  The  water  is 
filtered  by  Cole-Palmer  ion  X-change  cartridges  in  order  to  remove  any 
impurities  which  might  be  present  in  the  water.  The  water  pressure  in 
the  laboratory  is  sufficient  to  fill  the  tube  thus  rendering  the  pumping 
of  water  into  the  tube  unnecessary.  The  water  flows  readily  around  the 
transducer  plate. 

B.  Experimental  Procedure 

The  laboratory  sediment  was  prepared  by  adding  the  glass  beads  to  the 

water  rather  than  by  adding  water  to  the  glass  beads.  The  sediment  was 
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prepared  in  this  way  because  it  has  been  reported  that  glass  beads  have 
a  tendency  to  trap  large  amounts  of  gas  in  their  pores.  Experiments 
carried  out  on  small  samples  indicated  that  if  the  glass  beads  were  added 
to  the  water  no  significant  pockets  of  gas  were  trapped  within  the  body 
of  the  sediment.  When  water  was  added  to  the  glass  beads  large  pockets 

of  gas  were  formed  which  were  extremely  difficult  to  remove. 

The  sediment  was  prepared  by  introducing  a  quantity  of  water  into  the 
sound  tube  and  then  slowly  pouring  glass  beads  into  the  tube  from  the 

top.  The  walls  of  the  sound  tube  were  tapped  periodically  with  a  rubber 

mallet  to  help  ensure  an  even  packing  of  the  glass  beads. 

When  the  required  amount  of  sediment  had  been  prepared,  the  sound 
tube  was  sealed  by  placing  an  end-plate  at  the  top  of  the  tube.  The 


3.  The  electronic  equipment 


The  electronic  equipment  used  for  generating  and  detecting  tone 
bursts  is  shown  in  Fig.  1.  Sine  waves  were  produced  bv  a  Rockland  5100 
frequency  synthesizer.  The  Rockland  is  a  calibrated  frequency  source  so 
it  was  unnecessary  to  monitor  the  frequency  of  the  sine  waves  during  data 
acquisition.  Bursts  of  sine  waves  were  produced  by  passing  the  output  of 
the  Rockland  frequency  synthesizer  through  a  General  Radio  1396-B  tone- 
burst  generator.  The  tone-burst  generator  allowed  the  duration  of  the 
bursts  and  the  interval  between  the  bursts  to  be  controlled. 

In  the  experiments  the  duration  of  the  tone-bursts  was  maintained 
less  than  the  time  needed  for  the  burst  to  travel  twice  the  length  of  the 
tube  in  order  to  avoid  the  effect  of  standing  waves.  It  was  found  that 
it  was  impossible  to  set  up  standing  waves  in  the  tube  because  of  the 
high  attenuation  of  compressional  waves  in  sediments. 

The  tone-bursts  were  amplified  by  a  Dvnaco  Mark  VI  power  amplifier 
which  is  capacitively  coupled  to  the  dc  biased  transmitter.  The  signal 
received  by  the  hydrophone  is  amplified  by  a  Bruel  &  Kjaer  2626 
conditioning  amplifier.  The  output  of  the  conditioning  amplifier  is  then 
further  amplified  by  a  Princeton  Applied  Research  113  low-noise  pre¬ 
amplifier.  The  resulting  waveform  is  displayed  on  a  Nicolet  digital 
oscilloscope. 


SS T 


20  40  60  80 


■J  irm 


Figure  4 
Hydrophone  Plate 
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2.  The  receiver 

The  acoustic  signal  is  received  bv  a  Bruel  &  Kjaer  8103  miniature 
hydrophone.  The  B&K  8103  is  a  small,  high-sensitivitv  waterborne  trans¬ 
ducer  designed  for  making  absolute  sound  measurements  over  the  frequency 
range  0.1  Hz  to  200  kHz. 

The  hydrophone  is  mounted  in  a  threaded  brass  sleeve.  The  entire 
assembly  is  bolted  to  a  0.076  m  (3")  thick  stainless  steel  plate.  The 
hydrophone  plate  is  surrounded  by  a  Teflon  gasket  to  facilitate  its 
movement  in  the  sound  tube.  The  hydrophone  plate  sits  on  three  raising 
and  lowering  bolts  which  are  threaded  and  emerge  through  the  bottom  end- 
plate  of  the  sound  tube.  The  raising  and  lowering  bolts  allow  the  hydro¬ 
phone  plate  to  be  raised  or  lowered,  to  cover  or  expose  a  0.0254  m  (l") 
valve  through  which  the  sediment  can  be  pumped  out  of  the  sound  tube. 
The  hydrophone  assembly  is  diagrammed  in  Fig.  4. 

Before  the  hydrophone  plate  was  introduced  into  the  sound  tube  an 
experiment  was  performed  to  determine  the  minimum  distance  that  the 
hydrophone  could  protrude  from  the  front  of  the  hvdrophone  plate  without 
affecting  the  received  signal.  It  was  found  that  it  was  necessarv  for 
the  hydrophone  to  protrude  0.014  m  from  the  hvdrophone  plate;  this  means 
that  the  acoustic  center  of  the  hydrophone  is  0.004  m  above  the 
hydrophone  plate. 


caustic  solution  will  etch  aluminum  off  Mylar,  the  coating  of  paint 
prevented  this. 

The  transmitter  is  connected  to  a  second  aluminum  plate  by  means  of 

an  aluminum  spacer.  This  second  aluminum  plate  is  also  surrounded  by  a 

Teflon  gasket.  The  transmitter  when  assembled  has  been  described  bv 
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Shields  as  a  'spool-shaped  plunger. '  This  arrangement  reduces  the 
tendency  for  the  transmitter  to  twist  in  the  tube.  The  transmitter  is 
allowed  to  sink  to  the  level  of  the  sediment  under  its  own  weight.  This 
process  usually  takes  several  minutes  as  air  is  forced  around  the  edges 
of  the  Teflon  gasket.  The  transducer  is  raised  by  a  block  and  tackle 
suspended  from  the  top  of  the  scaffold. 

The  process  of  raising  and  lowering  the  transmitter  caused  most  of 
the  problems  encountered  in  acquiring  the  attenuation  data.  When  this 
research  project  was  originated  in  1979,  one  of  its  stated  aims  was  to 
adapt  expertise  in  determining  sound  attenuation  in  gases  to  laboratory 
sediments.  That  goal  has  been  achieved,  but  problems  became  evident  when 
the  large-scale  apparatus  was  used.  This  research  has  demonstrated  the 
viability  of  using  a  solid  dielectric  transducer  as  the  sound  source,  but 
the  experience  gained  in  using  the  system  calls  into  auestion  the  advisa¬ 
bility  of  using  a  moveable  sound  source.  The  transmitter  had  a  tendencv 
to  become  stuck  in  the  tube  because  the  glass  beads  tended  to  coat  the 
walls  of  the  sound  tube.  The  transmitter  was  not  very  durable,  several 
were  used  during  the  course  of  the  experiment. 


1.  The  sound  transmitter 


The  sound  transmitter  used  in  this  experiment  is  a  larger  version  of 
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the  solid  dielectric  capacitor  transducers  designed  by  Kuhl;  it  is  a 

waterproof  version  of  the  large  capacitor  microphones  used  previously  in 
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this  laboratory  by  Shields. 

The  transmitter  is  constructed  using  a  sheet  of  0.5-mil  Mylar  which 
is  aluminized  on  one  side.  The  Mylar  is  stretched  over  an  aluminum 
backing  plate  with  the  aluminized  surface  on  the  outside  electrically 
insulated  from  the  backing  plate.  A  dc  polarizing  voltage  is  applied 
between  the  aluminized  Mylar  and  the  backing  plate  and  an  ac  voltage  is 
applied  in  series  with  the  polarizing  voltage  by  means  of  the  voltage 
adder  circuit  shown  in  Fig.  2. 

The  transducer  was  designed  to  completely  fill  the  cross-sectional 
area  of  the  tube.  The  construction  details  of  the  transducer  are  shown 
in  Fig.  3.  The  Mylar  film  is  attached  to  a  floating  tension  ring,  and 
the  Mylar  is  stretched  over  the  backing  plate  by  means  of  three  tension 
screws.  The  backing  plate  is  coupled  to,  but  is  electrically  insulated 
from,  a  0.052  m  (2")  thick  stainless  steel  loading  plate.  The 
transmitter  is  surrounded  by  a  Teflon  gasket  which  was  designed  to  allow 
the  transmitter  to  fit  tightly  in  the  sound  tube  and  still  allow  it  to 
move  freely.  The  front  surface  of  the  aluminized  Mylar  was  coated  with 
an  alkvd  resin  based  paint;  the  oaint  water-sealed  the  transducer  and 
also  electrically  insulated  the  aluminized  Mylar  from  the  water.  When 
small  glass  beads  are  added  to  water,  they  form  a  caustic  solution;  a 


Block  Diagram  of  System 
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The  mean  fluid  pressure  in  the  tube  was  calculated  to  be  1.3  x  10 J  Pa. 

Unless  otherwise  noted  all  data  was  obtained  from  the  manufacturer's 
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specifications  for  the  gl-’ss  beads  or  Kinsler  and  Frey.  When  the  wall 

losses  were  calculated  for  the  sound  tube  containing  the  sediment  it  was 
found  that  the  viscous  losses  dominated  those  due  to  thermal  conduction 
at  the  walls. 

The  attenuation  in  the  sediment  with  the  wall  losses  subtracted  is 

presented  in  Fig.  6.  The  attenuation  measured  in  the  laboratory  sediment 

is  not  a  linear  function  of  frequency,  as  predicted  by  Hamilton;  ^  it 

varies  as  the  square  of  frequency  at  low  frequency  (below  4  kHz)  and 

shows  an  approximately  linear  dependence  on  frequency  above  4  kHz.  The 

Biot  model  predicts  that,  at  low  frequencies,  the  attenuation  of  compres- 

sional  waves  in  a  sediment  varies  as  the  square  of  frequencv,  and  at  high 

frequencies  the  attenuation  varies  as  the  square  root  of  frequency 

resulting  in  an  approximately  linear  variation  with  frequencv  at 
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intermediate  frequencies.  The  attenuation  measurements  reported  in 

this  dissertation  were  performed  at  low  and  intermediate  frequencies  for 
the  size  of  particles  used  in  the  laboratory  sediment.  As  explained  in 
Chapter  II  the  terms  low  and  high  frequency  are  dependent  on  the  size  of 
the  particles  in  the  porous  medium  when  the  Biot  model  is  used. 


Attenuation  (dB/m) 


Frequency  (kHz) 

Figure  6 

Attenuation  Corrected  for  Wall  Losses  Compared  to  the  Predictions 
of  the  Biot  and  Hamilton  Models 


VI.  CONCLUSIONS 


Although  the  functional  dependence  of  the  attenuation  measurements 
(precisely  in  accord  with  the  Biot  model)  is  gratifying,  the  actual 
magnitude  of  the  attenuation  above  4  kHz  is  higher  than  expected.  Figure 
6  contains  a  plot  of  Eq.  I,  the  attenuation  predicted  by  the  Hamilton 
model  using  the  maximum  value  of  the  parameter  K  (0.9)  for  particles  on 
the  order  of  65  microns  in  size.^  Hamilton's  model  predicts  an  attenua¬ 
tion  which  is  about  1  dB/m  lower  than  measured  here  in  the  range  from  5- 
10  kHz.  This  excess  attenuation  may  be  accounted  for  in  one  of  two  ways; 
the  method  of  calculating  the  attenuation  due  to  the  tube  walls  may  be 
inadequate  due  to  the  finite  rigidity  of  the  sediment,  or  the  sediment 
may  have  contained  some  gas  bubbles  the  presence  of  which  will  contribute 
to  the  measured  attenuation,  or  a  combination  of  both. 

The  relatively  high  attenuation  in  the  system  made  it  impossible  to 
determine  the  propagation  velocity  of  the  compressional  waves  using  the 
resonance  method.  The  propagation  velocity  is  useful  to  obtain  an  esti¬ 
mate  of  the  magnitudes  of  the  bulk  and  shear  moduli  of  the  frame.  The 
absence  of  this  estimate  permits  a  considerable  latitude  in  the  choice  of 
these  quantities  to  fit  the  experimentally  measured  attenuation  to  the 
Biot  model.  The  method  used  to  calculate  the  attenuation  predicted  bv 
the  Biot  model  was  to  choose  a  value  from  Poisson's  ratio  (which  relates 
the  shear  and  bulk  moduli  of  the  frame) and  then  to  varv  the  bulk 


modulus  of  the  frame  until  the  experimental  data  fits  the  theoretical 
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curve.  Fine  adjustments  are  made  by  altering  the  value  of  Poisson's 
ratio  when  the  experimental  data  lies  close  to  the  theoretical  curve.  In 
this  way  the  experimental  data  can  be  made  to  lie  as  close  as  desired  to 
the  theoretical  curve.  The  attenuation  predicted  by  the  Biot  model 
contained  in  Fig.  6  was  calculated  using  the  computer  program  at  Appendix 
A. 

The  functional  dependence  of  the  attenuation  data  displayed  in  Fig.  6 
is  a  more  convincing  testament  to  the  veracity  of  the  Biot  model  than  is 
the  theoretical  Biot  curve  in  the  absence  of  realistic  estimates  of  the 
frame  moduli. 

Future  attenuation  measurements  using  the  system  designed  for  this 

experient  will  be  facilitated  if  the  engineering  improvements  in  the 

transmitter  described  in  Chapter  IV  are  made.  Since  the  transmitter  was 

designed  and  constructed,  better  designs  for  underwater  solid  dielectric 
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transducers  have  appeared  in  the  literature.  Experience  gained  during 
the  course  of  this  investigation  suggests  that  it  would  be  advantageous 
to  fill  the  tube  with  water  above  the  level  of  the  sediment  and  to 
correct  the  resulting  measurements  for  the  (small)  attenuation  due  to 
water.  A  system  to  measure  the  propagation  velocity  of  the  compressional 
waves  would  enable  the  experimental  data  to  be  more  easily  fitted  to  the 
Biot  model,  but  it  is  recommended  that  methods  be  developed  to  directly 
determine  the  frame  moduli  to  eliminate  the  free  parameters  in  the  Biot 
model.  Future  theoretical  work  should  concentrate  on  either  the  calcula¬ 
tion  of  H,  C,  and  M  directly,  or  a  satisfactory  method  of  calculating  the 
frame  moduli  for  unconsolidated  media. 


In  conclusion  it  has  been  demonstrated  that  the  Biot  model 
successfully  predicts  the  attenuation  of  low-f recuencv  compressional 
waves  in  a  laboratory  sediment.  The  Biot  model  continues  to  pass  everv 
test  to  which  it  is  put  and  it  is  recommended  that  it  be  adopted  as  the 
standard  model  for  wave  propagation  through  ocean  sediments. 
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This  program  was  used  to  calculate  the  attenuation  of  sound  in  sedi¬ 
ments  using  the  Biot  model. 


C  Computation  of  sound  velocity  and  attenuation  in  sediments  using  the 
C  Biot  Model 
C  Reconstructed  version 
C 

C************* ********************************************* ************** 

C  TYPE  DECLARATIONS  OF  VARIABLES 

C************** *********************************************** *********** 

c 

REAL  KR,KF, PHI, PERM, VISC.NU.DIA 
REAL  RHO , RHOF , RHOS , RHOC 
REAL  PV1 ,PV2,ABS1,ABS2 
REAL  FREQ,AF, AF2 , AF3 , AF4 ,PI 
REAL  T0RT,X,QDC1 
COMPLEX  KB,MU,H,C,M 
COMPLEX  KRB , DKB , D , I , F 
COMPLEX  RPL.RMI ,Q1 ,Q2 
COMPLEX  QDA , QDB , QDB 1 , QDB 2 , QDC , QDC 2 
C 

I=(0. 0,1.0) 

1=3- 1415926 

READ  ( 1 0 , * ) KR , KF , KB ,MU , V I SC , PERM ,TORT , DI A , RHOS , RHOF , PHI 
Q  ********  k  kkkkkkkkkkkkkkkkkkkkkkk  k  *********************  *  ****************  * 

C  Calculation  of  H,C,M 

Qkkkkkkkkk  k  it  kit  irk  kirk  kirk  kit  it  irk  kirk  kick  it  kit  irk  irk  kit  irk  kick  kirk  kit  irk  irk  irk  kirk  kirk  kirk  irk  idc  irk 

c 

D=KR* ( 1 . 0+PH I * ( KR-KF ) /KF ) 

KRB=KR-KB 

DKB=D-KB 

H=(KRB*KRB/DKB)+KB+(4.0/3.0)*MU 

C=KR*KRB/DKB 

M=KR*KR/DKB 


C****************************************** ************** **************** 

C  Calculation  of  densities 

C k kkkkkkkkkkkkkkkkkkkk k k kkkkkkkkkkkkkkkkkk  kkkkkkkkkk  kkkkkkkkkkkk k k k kkkkk k 

c 

RHO=( 1.0-PHI )*RH0S+PHI*RH0F 
RH0C=T0RT*RH0F / PHI 
C 
C 
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C 

0* ******* ************ ******** ********* *********************************** 

C  Calculations  using  the  Biot  model 

q***** ************************************************************* ****** 

C 

QDA=(C*C)-(H*M) 

NU-VISC/PERM 

QDB 1 =M*RH0+H*RH0C-2 . 0*C*  RHOF 
QDC1=RH0F*RH0F-RH0*RH0C 
DO  1000  IML-1 ,28 
READ( 1 1 , *)FREQ 
AF»(PI*FREQ)+(PI*FREQ) 

AF2-AF*AF 

AF3-AF*AF2 

AF4«AF2*AF2 

X-AF*RHOF/VISC 

X=SQRT(X)*DIA/7.0 

CALL  DYNVIS(X.F) 

QDB2-I*H*NU*F 
QDC2«I*RHO*NU*F 
QDB=(QDB1*AF2 )-(QDB2*AF) 

QDC»(QDC1*AF4)+(QDC2*AF3) 

CALL  SOLVE (QDA, QDB, QDC,RPL,RMI) 

RPL-CSQRT(RPL) 

RMI-CSQRT(RMI) 

0 ******** ****** ******************************************* *************** 

C  Calculation  of  phase  velocity  and  attenuation  in  nepers/m 

0* ******** ******** ************************************************* ****** 
P V 1 =AF / RE AL ( RMI ) 

PV2=AF/REAL(RPL) 

ABS 1=-IMAG(RMI) 

ABS2— IMAG(RPL) 

WRITE ( 12,*)FREQ,PVl ,ABS1 
WRITE ( 1 3 ,*)FREQ,PV2 , ABS2 

1000  CONTINUE 
END 

Q ******** ****************************************************** ********** 

C  SUBROUTINE  DYNVIS(X,F) 

C  Routine  to  calculate  Biot's  dynamic  viscosity  factor 

0****** ********* ********* ***************************************** ******* 
SUBROUTINE  DYNVIS(X.F) 

COMPLEX  I ,T ,F 

IF  (X.LE.0.1)  GO  TO  2000 

BER-1.0 

BERP-0.0 

BEI=0. 0 

BE IP-0.0 

I-(0. 0,1.0) 

C 


DO  1001  K=1 , 5 
N-K 

A=(-1.0)**N 
BRE2=2.0*N 
BRE4=4. 0*N 
BBER-X**BRE4 
CBRE»2.0*BRE4 
DBRE=GAMMA( BRE2+1 . 0) 

EBRE-DBRE*DBRE 

TBER”( A*BBER) / (CBRE*EBRE ) 

BER-BER+TBER 

BERP4M«4.0*N-1.0 

BBERP-X**BERP4M 

TBERP-(A*BBERP*BRE4)/(CBRE*EBRE) 

BERP-BERP+TBERP 

1001  CONTINUE 

DO  1002  K-1,6 

N-K-l 

A— 1.0**N 

BEI2-2.0*N+1.0 

BEI4»4.0*N+2.0 

BBEI»X**BEI4 

CBEII»2.0*BEI4 

DBEII«GAMMA(BEI2+1.0) 

EBEII«DBEII*DBEII 

TBEI-(A*BBEI)/(CBEII*EBEII) 

BEI-BEI+TBEI 

BEIP4M=4.0*N+1.0 

BBEIP»X**BEIP4M 

TBEIP»(A*BBEIP*BEI4)/(CBEII*EBEII) 

BEIP-BEIP+TBEIP 

1002  CONTINUE 

C******************************* *************************** ************** 

C  Computation  of  dynamic  viscosity  factor 

C*** ************ ****************************************** *************** 

T=(BERP+I*BEIP)/(BER+I*BEI) 

F=( 0 . 2  5*X*T) / ( 1 . 0-( ( 2 . 0*T) / ( I*X ) ) ) 

RETURN 

2000  F«(l. 0,0.0) 

RETURN 


END 

SUBROUTINE  SOLVE ( A , B , C , RPL , RMI ) 
COMPLEX  A,B,C,RPL,RMI,DISCR 
DISCR=CSQRT ( B*B-4 . 0*A*C ) 
RPL=(-B+DISCR)/(2.0*A) 
RMI-(-B-DISCR)/(2.0*A) 

RETURN 


END 
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